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ABSTRACT 

We study the stability of charged dust grains orbiting a planet and subject to gravity and the 
electromagnetic force. Our numerical models cover a broad range of launch distances from the 
planetary surface to beyond synchronous orbit, and the full range of charge-to-mass ratios from 
ions to rocks. Treating the spinning planetary magnetic field as an aligned dipole, we map regions 
of radial and vertical instability where dust grains are driven to escape or crash into the planet. 
We derive the boundaries between stable and unstable trajectories analytically, and apply our 
models to Jupiter, Saturn and the Earth, whose magnetic fields are reasonably well represented 
by aligned dipoles. 
Introduction 



The discoveries of the faint dusty ring systems 
of the giant planets beginning in the late 1970s 
greatly changed our understanding of planetary 
rings. Unlike Satu rn's classical rings, which are 
most likely ancient (|Canudl2010l) . dusty rings are 
young and are continually replenished from source 
satellites. Individual ring particles have short life- 
times against drag forces and other loss mecha- 
nisms, and because dusty rings are so diffuse, they 
are essentially collisionless. Furthermore, dusty 
rings are affected by a host of non-gravitational 
forces including solar radiation pressure and elec- 
tromagnetism, which can sculpt them in interest- 
ing ways. 

Since the giant planets are far from the Sun 
and dusty rings are normally near their primary, 
radiation pressure is usually a weak perturbation 
to the planet's gravity. The electromagnetic force 
arising from the motion of charged dust grains rel- 
ative to the planetary magnetic field, however, can 
be quite strong. In particular, with nominal elec- 
tric charges, dust grains smaller than a fraction of 
a micron in radius are more strongly affected by 
electromagnetism than gravity. 

Dust in space acquires electric charges in several 
ways. Moving through the plasma environment 
produces a negative charge on a grain, since the 
plasma electrons are much lighter and swifter than 



ions and hence are captured mor e frequently by 
orbiting dust grains ([Goertz 1989). On the other 
hand, sunlight ejects photo-electrons from the sur- 
face of a grain, and can cause positive charges 
(jHoranvi et Electron or ion impacts will 

also produce secondary electron emission, which 
also favor s a net positiv e equilibrium charge on 
the grain (|Whipplelll98l . These currents interact 
in complicated ways; the charging of a dust grain 
depends on the physical properties of the grain it- 
self a n d also on its charg e history (jMever-Vernet 
1982 ) . Graps et all ( 20081 ) provide an excellent re- 



view of these processes. 

Many authors have investigated detailed as- 
pects of the motion of charged grains in planetary 
magnetic fields, but no study has yet determined 
the orbital stability of grains for all charge-to-mass 
ratios launched at all distances in a systematic 
way. In this paper we explore the local and global 
stability of both positive and negative dust grains 
launched from ring particle parent bodies which 
themselves orbit at the local Kepler speed. 

1.1. Motion in the Kepler and Lorentz 
Limits 

As grains with radii greater than several mi- 
crons have small charge-to-mass ratios, electro- 
magnetic effects are weak, and the grains orbit 
the planet along nearly Keplerian ellipses. In the 
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frame rotating with the mean motion of the dust 
particle, the orbits appear as retr ograde ellipti 



cal ep icycles with a 2:1 aspect ratio (|Mendis et al. 
2J). When gravity acts alone, the vertical, ra- 



ror p oi nts north and south of th e equa tor (jStormer 
19551) . iThomsen and van Alien! (|l980f ) studied the 
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dial and azimuthal motions all have precisely the 
same frequency. Equations governing the slow 
changes to the ellipse's orbital elements due to 
weak electromagnetic perturbations from a rotat- 
ing aligned dipo le magnetic field are given by 
Hamilton! (|l993al ). These equations show that the 
three frequencies diverge slightly and are functions 
of the sign and magnitude of the charge as well as 
the distance from the planet and synchronous or- 
bit. 

Conversely, the very smallest dust grains ap- 
proach the Lorentz limit, where the electromag- 
netic force dominates over gravity. In this regime, 
the frequencies of radial, vertical and azimuthal 
motions differ significantly. The radial oscillation 
is fastest and, as the electromagnetic force is per- 
pendicular to the rotating magnetic field, particles 
gyrate about local field lines on typical timescales 
of seconds for dust, and microseconds for ions. 

Dust grains typically oscillate vertically on a 
timescale of hours to days. Since this timescale 
is far slower than gyration, an adiabatic invariant 
exists and can easily be found. In the absence of 
forces other than electromagnetism, and ignoring 
planetary rotation, a dust grain's speed v remains 
constant: v' 2 = v\ + ujj , where v± and «|| are the 
speeds perpendicular and parallel to the magnetic 
field lines, respectively. The v±_ component de- 
termines the radius of the gyrocycle, while the v\i 
component moves the center of gyration to regions 
of differing magnetic field strength. If changes 
to a non-rotating magnetic field B are small over 
the size and time scales of gyromotion, the ratio 
Vj_/B, where B is the loca l strength of the field, 



is an adiabatic invariant (jde Pater and Lissauer 
201dh and hence is nearly constant. These two 
conditions provide an important constraint on the 
grain's motion parallel to the field lines. As a 
grain with a vertical velocity component climbs 
up a magnetic field line away from the equato- 
rial plane, the field strength B increases, v± also 
increases, and hence v\\ must decrease. There is 
thus a restoring force towards the equatorial plane 
where the magnetic field strength is a local min- 
imum, and the motion parallel to the field lines 
takes the form of bounce oscillations between mir- 



bounce motion of particles in the Lorentz limit at 
Saturn. Their results neglected the effects of plan- 
etary rotation, and hence are most applicable to 
slow rotators like Mercury and potentially some 
planetary satellites. 

Finally, on the longest timescales (days), parti- 
cles drift longitudin ally with respect to the rotat - 
ing magnetic field ( de Pater and Lissauer! [201oh . 
forced by a number of effects including gravity, 
the curvature of the magnetic field, and VB. Be- 
cause these motions are usually slow compared to 
the gyration and bounce frequencies, it is often 
useful to assume that in the Lorentz limit, grains 
are tied to the local field lines. 

1.2. Dust Affected by both Gravity and 
Electromagnetism 

For a broad range of grain sizes from nanome- 
ters to microns, both gravity and the Lorentz 
force are significant, and their combined effect 
causes a number of dynamical phenomena that 
are distinct from either limiting case. As dust 
in this siz e range predominates in many plane- 
tary rings (iBurns et qJl999Hde Pater et all\ 19991 



IShowalter et all 120081 : kriiger et all l2009h . their 
dynamics have attracted much attention. 

ISchaffer and Burns! l|l994h provide a general 
framework for the motion of dust started on ini- 
tially Keplerian orbits. Since the radial forces 
on a dust grain at launch are not balanced as 
they are for a large parent body on a circular 
orbit, these dust grains necessarily have non- 
zero amplitude epicyclic motion. For the mag- 
netic field configurations of the giant planets, 
a negatively-charged dust grain gyrates towards 
synchronous orbit while positively-charged dust 
initially moves away from this location. In fact, 
some positively-charged grains are radially unsta- 
ble and either crash into the planet if launched 
inside synchronous orbit, or are expelled outwards 
if launched from beyond this distance. The lat- 
ter have been detected as high-speed du st streams 
near Jupiter dGriin et all 1 199a Il998l) and Sat- 



urn (jKempf et all I2005T ). Theoretical explana- 
tions for the electrom agnetic acceler a tion p ro- 
cess have been given by Horanvi et a l l dl993al!3), 
lHamilton and Burns! ( 1993b! ) and Graps et all 
(!2000h . 
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Mendis et al. I (GUI) explored the shape and fre- 
quency of epicycles for negatively-charged grains 
in the transitional regime, where both EM effects 
and gravity are comparable. The epicycles make 
a smooth transition from perfectly circular clock- 
wise (retrograde) gyromotion in the Lorentz limit, 
to 2:1 retrograde elli p tical epicycles in the Kepler 
limit. iMitchell et al. ( 2003 ') studied the shapes of 
epicyclic motion for positive grains and found that 
there is not a similarly smooth transition from pro- 
grade gyromotion to retrograde Kepler epicycles, 
and that the epicyclic motions of intermediately- 
sized grains cannot be represented as ellipses. The 
effects of gravity and electromagnetism compete 
for intermediate charge-to-mass values and motion 
can be primarily radial, leading to escape or col- 
lision ( Horanvi et al. 1993a ; Hamilton and Burns! 
1993b!) . 

iNorthrop and Hilll (|l982l Il983ah and Northrop 
and Connerney (1987) studied the vertical mo- 
tion of negatively-charged dust grains on circu- 
lar uninclined orbits in a centered and aligned 
dipole field, a configuration most closely realized 
by Saturn. They found that some small grains 
on initially centrifugally-balanced circular trajec- 
tories inside the synchronous orbital distance are 
locally unstable to vertical perturbations, climb- 
ing magnetic field lines to crash into the planet 
at high latitudes. Some motions at high lati- 
tude, however, are stable: iHoward et al. ( 19991 
2000h identified non-equatorial equilibrium points 
for charged dust grains, and showed than dust 
grains can orbit them stably. They character- 
ized these "halo" orbits for positive and nega- 
tive charged grains on both prograde and retro- 



grade trajectories. Howard and Horanvil (|200ll ) 



used these analytical results to argue for a sta- 
ble population of positively-charged grains in ret- 
rograde orbits and developed numerical models of 
such halo dust populations at Saturn. Grains that 
may populate these halos, however, are unlikely 
to result from the equatorial launches considered 
here. 

If one of the dust grain's natural frequencies 
matches a characteristic spatial frequency of the 
rotating multipolar magnetic fi eld, the particle ex- 
periences a Lorentz resonance dBurns et aLlll985t 
Schaf fei~and Burn"s1ll987i Il992t lHamilton and Burns 



and can have a dramatic effect on a dust grain's or- 
bit, exciting large radial and/or vertical motions. 
These resonances have been primarily studied in 
the Kepler limit appropriate for the micron-sized 
particles seen in the dusty rings of Jupiter. In our 
idealized problem, with an axisymmetric magnetic 
dipole, Lorentz resonances cannot occur. 

Variations in a dust grain's charge can also alter 
its trajectory over surprisingly rapid timescales. 
Gradients in the plasma properties, including den- 
sity, temperature and even composition affect the 
equilibrium potential of a grain by altering the 
direct electron and ion currents. This can re- 
sult in resonant charge variation with gyrophase, 
causi ng radial drift. Working in the Lorentz 
limit, INorthrop and Hilll (|1983al) noted that with 
large radial excursions, the grain's speed through 
the plasma can vary significantly with gyrophase, 
leading to enhanced charging at one extremity. 
A similar effect occurs in the Kepler limit where 
resonant charge variation can cause a dramatic 
evolution in the orbital e l ements of a dust grain 
(|Burns and Schafferll 19891 ). INorthrop et q/.l(|l989h 
found that the varying charge has a time lag 
that depends on the plasma density and grain 
capacitance. These time lags can cause grains 
to drift towards or away from synchronous or- 
bit depending on the grain speed, and on any 
radial te mperature or dens it y gra dients in the 
plasma. ISchaffer and Burns! (|1995l ) explored the 
effects of stochastic charging on extremely small 
grains, where the discrete nature of charge cannot 
be ignored. They found that Lorentz resonances 
are robust enough to survive even for small dust 
grains with only a few electric charges. 

The dynamics of time-variable charging may 
play an important role in determining the struc- 

ture of Saturn's E ring ( Juhasz and Horanvill2004 ) 

and J upiter's main ring and halo (jHoranvi and Juhasz 
2010h . Another example of charge variation occurs 



when the insolation of a dust grain is interrupted 
during transit through the planetary shadow. This 
induces a variation in charge that resonates with 
the g rai n's orbital frequency ( Horanvi and Burnsl 



199lh . lHamilton and Kriigerl (|2008l ) found that 



1993a! : lHamilton! 1 1994) . Lorentz resonances be- 
have similarly to their gravitational counterparts 



this shadow resonance excites radial motions while 
normally leaving vertical structure unaltered. 
This effect can explain the appearance of the faint 
outward extension of Jupiter's Thebe ring, and 
the properties of its dust population sampled by 



3 



the Galileo dust detector ( Kriiger et al. 2009) 



1.3. Research Goals 

In this study, we consider the orbits of charged 
grains launched in planetary ring systems. Our 
aim is to explore the boundaries between stable 
and unstable orbits in aligned and centered dipolar 
magnetic fields. Dipolar fields have the advantage 
of being analytically tractable while still captur- 
ing most of the important physics. Under what 
conditions are grains unstable to vertical pertur- 
bations? Which grains escape the planet as high 
speed dust streams? And which grains will strike 
the planet after launch? All of these instabilities 
depend on the launch distance of the grain and 
its charge-to-mass ratio. We first explore grain 
trajectories numerically and then derive analyti- 
cal solutions for the stability boundaries that we 
find. 

There are several standard choices for express- 
ing the ratio of the Lorentz and gravitational 
forces. The charge-to-m ass ratio q/m in C/kg 
( Northrop and Hill Il982j) or in statCoulomb/g 
" 20031) ma 



([Mitchell et aZ.1 120031) may be the most straight- 
forward, but it is cumbersome. For this reason, 
converting to the grain potential measured in 
Volts, which is constant for different-sized dus t 
is a common choi ce ( Mendis et all 198S 
Howard et all I200C 



grains 
Schaffer and Burns 



1994; 



Mitchell et all 120031 ). Yet another option is to 



express the charge-to-mass ratio in terms of fre- 
quencies associated with the primary motions of 
the grain, such as the gyrofrequency, orbital fre- 
quency and the sp i n frequency of the p lanet (eg. 
Mendis et aZlll982t iMitchell et aLll2003f) . 

We choose a related path, namely to fold q/m 
and key planetary parameters i nto a single dimen - 
sionless parameter L* following iHamiltonl ( 1993a| ). 
Consider the Lorentz force in a rotating magnetic 
field: 

F B = -{v - fl x r) x B, (1) 
c 

where c is the speed of light, r and v are the grain's 
position and velocity in the inertial frame, O is 
the spin vector of the planet, and B is the mag- 
netic field. We use CGS units here and through- 
out to simplify the appearance of the electromag- 
netic equations. The second component of Eq. [T] 
is qE, where E = — -(Q x r) x B is the so-called 
co-rotational electric field which acts to accelerate 



charged grains across magnetic field lines. Since a 
dipolar magnetic field obeys B = ~gioRp/r 3 z in 
the midplane (with gio the magnetic field strength 
at the planet's equator), E, like gravity, is propor- 
tional to 1/r 2 there. Thus the ratio of the electric 
force to gravity is both independent of distance 
and dimensionless: 



qgwRptt 

GMpTTlC 



(2) 



Here, R p and M p are the planetary radius and 
mass, m is the dust grain mass and G is the grav- 
itational constant. Note that the the sign of L* 
depends on the product of two signed quantities, 
q and gio. For all of the giant planets, the mag- 
netic north pole is in the northern hermisphere, 
and gio > 0. However, for the Earth at the current 
epoch, 510 < and the magnetic and geographic 
poles are in opposite hemispheres. 

We have made a s light n otational change 
L — > from Hamilton (|l993alibh to avoid confu- 
sion with the L-shell of magnetospheric physics. 
Choosing as an independent variable takes the 
place of assuming a particular electric potential, 
grain size and grain density. We focus our study 
primarily on Jupiter, the planet with by far the 
strongest magnetic field, but also apply our results 
to Saturn and to the Earth. 

2. Numerical Simulations 

Approximating Jupiter's magnetic field as an 
aligned dipole by including just g w = 4.218 Gauss 
(|Desslerlll983l ). we tested the stability of dust grain 
orbits over a range of grain sizes and launch dis- 
tances both inside and outside synchronous or- 
bit. We used a Runge-Kutta fourth-order inte- 
grator and launched grains at the local Kepler 
speed with a small initial latitude of A = 0.01°. 
This tiny nominal value ensures a launch close to 
the midplane, whilst avoiding potential numerical 
problems of launching a grain precisely at A = 0. 
Non-zero launch speeds from the parent particle 
do have a small effect on the stability boundaries, 
one that we will explore in more depth in a future 
study. 

Our models treat the grain charge as constant 
and neglect J2, other higher-order components of 
the gravitational field, and radiation pressure. For 
both negative and positive grains, we ran sim- 



4 



ulations for a grid of 80 values of and 100 
launch distances (tl)- The charge-to-mass ratio 
spans four decades from the Lorentz regime where 
EM dominates (|L*| >> 1), to the Kepler regime 
where gravity reigns (|L*| << 1). The range of 
launch distances extends from the planetary sur- 
face to well beyond the synchronous orbital dis- 
tance (R syn ), and trajectories were followed for 
up to 0.1 years. With some experimentation, we 
determined that all relevant dynamical timescales 
are < 0.1 years and that for longer integration 
times, the appearance of our stability plots does 
not change significantly. 

In Fig. [1] we plot the fate of 8000 negative and 
8000 positive dust grains and find complex regions 
of instability. The negatively-charged dust grains 
in Fig. QJi display only vertical instability at mod- 
erate to high and inside R syn - Some are bound 
by high latitude restoring forces (locally unstable, 
light grey) whilst others crash into the planet at 
high latitude (both locally and globally unstable, 
darker grey). To separate these globally stable 
grains from locally stable ones, we choose a lati- 
tude threshold at A TO = 5°. Although 5° is a small 
latitude, it is far greater than the launch latitude 
of 0.01°; any grains excited beyond X m are clearly 
locally unstable, and we determined that our re- 
sults were fairly insensitive to actual value of A m . 



a d (|im) 



Northrop and Hilll (|1982l ) derived a boundary 



for the threshold between locally stable and un- 
stable trajectories for negatively-charged dust and 
found that grains launched within a certain dis- 
tance should leave the equatorial plane (NH82 
curve in Fig. [TJl) . In the Lorentz limit, the ver- 
tical instability allows grains to climb up local 
magnetic field lines into regions of stronger mag- 
netic field, while for smaller L* the path taken by 
these grains follows the lines of a pseudo-magnetic 
field w hich includes the effects of planetary ro- 
tation (jNorthrop and Hilll [1983) . The Northrop 
curve however, is not a good match to our data 
which reveal additional stable orbits (white areas) 
immediately inside this boundary and also close to 
the planetary surface. These diffe r ences arise from 
the fact that iNorthrop and Hilll (1982) assumed 
that grains are launched at their equilibrium cir- 
cular speeds, which differ from the circular speeds 
of parent bodies when ^ 0. Conversely, we 
launch our grains at v = J GM p /r, the circular 
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Fig. 1. — Stability of Kepler-launched a) (neg- 
ative) and b) (positive) dust grains at Jupiter. 
We model the planet with a spherically-symmetric 
gravitational field, and a centered and aligned 
dipolar magnetic field. All grains were launched 
with an initial latitude of A = 0.01° and followed 
for 0.1 years. The horizontal dashed line in both 
panels denotes the synchronous orbital distance 
at R syn — 2.2AR p . The grain radii (ad) in mi- 
crons along the upper axis are calculated assuming 
a density of 1 g/cm 3 and an electric potential of 
±5V so that |L*| = 0.0284/^. Dust grains in the 
white regions and lightest grey areas survive the 
full 0.1 years, with the latter reaching latitudes A 
in excess of 5°. Grains in the moderately-grey ar- 
eas are vertically unstable and strike the planet, 
also at high latitudes (A > 5°). The darkest re- 
gions, seen only in panel b, are radially unstable 
grains that crash into the planet (those with rj, < 
Rsyn), or escape to beyond r esc = 30i? p (from 
ri > R syn ) at latitudes less than 5°. We overplot 
three anal ytically-derived stab i lity b oundaries, ob- 
tained by INorthrop and Hilll (Tl 982) for negative 
grains, by Horanvi et all (Il993al) for sm a ll posi - 
tive grains, and bv [Hamilton and Burnsl (|l993bh 
for large positive grains. Each point on the plot 
is a trajectory, some of which (marked by filled 
squares), are illustrated in detail in Figs. [2] to O 



speed of the parent body, which is appropriate for 
debris produced by cratering impacts into these 
objects. In section 5, we develop a vertical stabil- 
ity criterion appropriate for our launch conditions. 

The situation for positive grains is quite differ- 
ent. Figure [TJd shows a less extensive region of 
vertical instability than Fig. QJ,, and one that is 
not active close to Jupiter. More dramatic, how- 
ever, are two regions of radial instability (darkest 
grey areas), separated by the synchronous orbital 
distance. Grains inside R syn are driven to strike 
Jupiter, while those outside escape the planet. If 
grains move beyond r esc = 30i? p , the inner mag- 
netosphere, we consider them to have escaped. As 
with A m , our numerical results are fairly insensi- 
tive to the exact value chosen for r esc , so long as 
it is large. 

To characterize the individual trajectories that 
make up Fig. [TJ we explore a few examples in de- 
tail, focusing on the positively-charged dust grains 
and proceeding from smaller to larger grains. Fig- 
ure [5] shows the trajectory of a dust grain that 
becomes vertically unstable and crashes into the 
planet at high latitude. These smallest grains 
spiral up magnetic field li nes, which for a dipole 
are g iven by rj cos 2 A = (ide Pater and Lissauer 
2010); collision with the planet or reflection from a 
high latitude mirror point typically occurs within 
a few tens of hours. By contrast, Fig. [3] shows an 
electromagnetically-dominated grain that remains 
stable at low latitude. 

A more subtle interplay between radial and ver- 
tical motions is illustrated in Fig. 3J This grain is 
outside the radial instability region in which grains 
collide with the planet at low latitude (darkest 
grey). Instead, large radial motions lead to insta- 
bility in the vertical direction, and ultimately, the 
grain strikes the planet at high latitude. Notice 
the two white dots near (L* = 1.34, r^/Rp = 2.2) 
in Fig. [TJd, signifying grains that survive the full 
0.1 year integration. These trajectories are indeed 
stable (for at least 100 years) and, as the effect is 
much more prominent for the Earth, we discuss it 
in more detail in section 6. 

Finally, Fig. [5] shows a dust grain just inside 
the lHamilton and Burns! (|1993bf ) L* = ~ stability 
limit. Although the dust grain does not escape, 
the non-linearity of its radial oscillation is large 
enough to excite substantial vertical motions. 
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Fig. 2. — The trajectory of a positively-charged 
grain orbiting Jupiter after launch at tl = 1.74i? p , 
with L* = 31.31 (a d = 0.03^m). Wc plot the 
scaled distance and latitude of the dust grain 
against time. The small, rapid radial gyration is 
just visible in the upper plot. The dust grain is 
vertically unstable on a much longer timescale and 
ultimately crashes into the planet. This trajectory 
is the left-most filled square in Fig. [Tp. 
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Fig. 3. — The trajectory of a stable positively- 
charged grain orbiting Jupiter after launch at Tl = 
1.74Rp, with L* = 3.04 (a d = 0.097^m). The 
grain undergoes radial oscillations much larger 
than in Fig. [5] but its latitude remains low. Here 
the bounce period is ~ 7 times longer than the 
gyroperiod. 
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Fig. 4. — The trajectory of a positive grain inside 
R S yn [r L = 2.0i? p ,L„ = 1.908, a d = 0.122/xm). 
Here, unlike Fig. 02 large radial motions ultimately 
excite vertical motions, forcing the trajectory to 
end with a collision at the planetary surface after 
just a few orbits. 



A glance at Fig. [Q shows th at most stability 
bound aries are unexplained. The lNorthrop and Hill 
(|1982n vertical stability boundary does not match 
the numerical data especially well, and only ap- 
plies to nega ti ve gra ins. For positive grains, 
Horanvi et al. I (|l993bh provided an approximate 
criterion for radial escape, which they applied 
far from synchronous orbit near Io. Their crite- 
rion is based on a comparison between the ra- 
dius of gyromotion r g , and the length scale over 
which the magnetic field changes substantially, 
namely where \B/(r g VB)\ ps 10, with the gyro- 
radius calculated in the Lorentz limit. Although 
not intended for use near synchronous orbit where 
r g — > 0, we neverthe less plot it on the lef t side o f 
Fig.[lp. Finally, the lHamilton and Burnsl (jl993bl ) 
L* — \ limit, derived from an energy argument, is 
a good match to the largest escaping grains. There 
is however, no analytical model for the broad class 
of grains that strike the planet. Accordingly, we 
seek to develop a unified theory that can cleanly 
determine all of these boundaries. We take up this 
task first for radial and then for vertical motions. 
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Fig. 5. — The trajectory of a positive grain outside 
Rsyn (r L = 2.7i? p ,L» = 0.419, a d = 0.26^m). 
As in Fig. 0] large radial oscillations eventually 
excite large vertical oscillations. Since the dust 
grain has L* < it is energetically requ ired to 
remain bound (jHamilton and Burnslll993bi ). Here 
T is measured in Earth days. 



3. Local Radial Stability Analysis 

Consider a centered magnetic dipole field that 
rotates with frequency f2 aro und a vertical axis 
aligne d in the z-direction. iNorthrop and Hill 
(|1982l ) derived the Hamiltonian for a charged dust 
grain in the rotating frame in cylindrical coordi- 
nates: 



H = U(p,z) 



p 2 + z 2 



(3) 



where p and i are the radial and vertical velocity 
components. The potential is given by 



U(p,z) = 



1 

v 



_ GM p p 2 L, 

m 



GAL 



nr 3 

L^p 2 



- 1 



(4) 



whe re the spherical radius r satisfies r 2 = p 2 + 
z 2 (INorthrop and Hilll Il982t ISchaffer and Burnsl 



1994 iHoward et all \200& iMitchell et all 120031 ). 
Equation |4] is the sum of two energetic compo- 
nents: first the azimuthal specific kinetic energy, 
which can be expressed as a function of r using 
the conservation of angular momentum, and then 
the potential associated with both the corotational 
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electric field and gravity. Note that we have cho- 
sen the zero of our potential to be approached as 
p — >• oo. Because t/(p, z) is independent of <j), the 
azimuthal coordinate, the canonical conjugate mo- 
mentum is a constant of the motion. For our 
launch condition from a large parent body on a 
circular orbit at r = r^: 



m 



r 2 L (n L + n gL ) 



(5) 



(jSchaffer and Burndll994 . where Ul and fl g L are 
the Kepler frequency and gyrofrequency evaluated 
at the launch distance r^: 



and 



I GM 



qB 



- L* 



(6) 



(7) 



mc ft 

Notice that in the gravity limit (L» — > 0), Eq. [5] 
reduces to r^riL, the specific angular momentum 
about the planet, while in the Lorentz limit (£* — > 
oo), it is r^QgL, the specific angular momentum 
about the center of gyromotion that moves with 
the magnetic field. 

If the motion of the particle is radially stable, 
it exhibits epicyclic motion about an equilibrium 
point determined from Eq. HI The existence of 
equilibrium points requires th at 7^7 = 7/7 = 0, 



both in the equatorial plane dNorthrop and Hill 
1982 ) and at high latitudes ( Howard et al. 1999, 
2000). The local stability of the equilibrium 
points, defined as whether oscillations about 
these points remain small, is then determined by 
considering the second derivatives of the poten- 
tial. Given our launch condition, we focus on 
the equatorial equilibrium points which are of 

d 2 u 



greatest interest. For these, 



= 0, 



dpdz 

p=p c ,z—0 

r — > p, and radial and vertical motions are ini- 
tially decoupled and may be considered separa tely 
(jNorthrop and Hil]|[l98a iMitchell et aLll2003f) . 

The equilibrium point is the guiding center of 
epicyclic motion. Grains launched at the guiding 
center have canonical conjugate momenta that 
are different from our Kepler-launched grains: 



namely, 



P c {lo c + ilgc), where w c is the orbital 



frequency of a grain at the guiding center, f2 ffc is 
the gyrofrequency at the guiding center, and p c is 



the guiding center distance in the equatorial plane. 
A local radial stability analysis is most relevant 
for our Kepler-launched grains if an equilibrium 
point is not too distant. Accordingly, it is impor- 
tant to distinguish between quantities evaluated at 
the Kepler launch position and those determined 
at the guiding center. Here and throughout, we 
use the subscript c for the guiding center and the 
subscript L for the launch position. At the equi- 



librium point, ^ 



0, which evaluates 



to: 



u 2 c Pc 



1 — - 



GAL 



0. (8) 



Physically, Eq.|8]just implies a balance offerees in 
the rotating frame, whereby the centrifugal force, 
the Lorentz force and gravity sum to zero. We 
solve Eq. [8] for the angular speed of the guiding 
center uj c , and find two real roots for L* < 1, 
which includes all negative charges. For L* > 1 
conversely, two equilibrium points exist only if 



R 3 



< 



4(L* - 1) 



(9) 



Two equilibria always exist inside R syn and every- 
where for >> 1 and £* << 1. There are no 
equilibrium points in a region starting at (L* = 2, 
Pc — Rsyn) in Fig. \Tjp, and opening upward to 
include an increasing range of L* values for in- 
creasing distance p c . In this region, no equilibrium 
point exists, and grains are guaranteed to be lo- 
cally unstable. Not surprisingly, this region is fully 
contained within the unstable portion of Fig. [TJd 
(darkest grey region outside R syn ). The existence 
of an equilibrium point, therefore, is a necessary 
prerequisite for stability. 

Additional instability in Fig. [TJd comes from two 
sources: i) the intrinsic instability of the equilib- 
rium point, if it exists, and ii) large amplitude 
motions about a locally stable equilibrium point. 
Large oscillations are beyond the scope of a local 
stability analysis and so we focus on small am- 
plitude radial motion near an equilibrium point, 
which takes the form 



d 2 U 
dp 2 1 



0. 



(10) 



Small radial motions are stable when §^ 



p—p c ,z—0 



s 



Kg > 0, which, from Eq. [4] can be written as: 



with dimensionless coefficients 



K7. = tjJ 2 , — 4cj r £l„ 



(ii) 



Il982t Northrop and Hill 1982; 
2003). Note here that the gyrofre- 



(jMendis et all 
Mitchell et al. 
quency fl gc is evaluated at the guiding center, 
and is given by Eq. [7] with the subscript change: 
L — > c. The epicyclic frequency k c reduces to the 
Kepler orbital frequency n c at the guiding center 
r c in the gravity limit, and to the gyrofrequency 
fl gc in the Lorentz limit. Radial excursions in 
both of these cases are small and, since k% > 0, 
are guaranteed to be stable. 

Radial motions are also initially small near 
synchronous orbit where electromagnetic forces 
are very weak (Eq. [T]), and so a local stability 
analysis is also applicable. At synchronous or- 
bit, uj c — n c — and Eq. [TT] reduces to = 
ri 2 (1 — 4L* + Li), which is positive for small or 
large L*. For 2 — \/3 < L* < 2 + -s/3, however, 
Eq.[TU]shows that radial motions near synchronous 
orbit are locally unstable. Comparing this analysis 
with Fig.QJ), we see that all orbits with rj, ~ R syn 
that are locally stable are, not surprisingly, also 
globally stable. The converse, however, does not 
hold: although most of the locally unstable or- 
bits are also globally unstable, some are in fact 
globally stable {e.g. L* < | just outside R syn in 
Fig. QJ>). In conclusion, the local analysis is con- 
sistent with our numerical experiments but cannot 
fully account for our stability boundaries. Accord- 
ingly, we turn to a global analysis, pausing first to 
put the potential of Eq. [4] into a more useful form 
and to derive the radius of gyration, r g . 

3.1. Radius of Gyration 

With our launch condition, grains are often far 
enough from an equilibrium point that the small 
oscillation approximation of Eq. [I0]is invalid. This 
is particularly true far from R syn and for L, «1. 
Returning to the effective potential of Eq. @] with 
the canonical conjugate momentum determined by 
launching the grain at the Kepler speed (Eq. [SJ , 
and limiting our attention to planar orbits for 
which z — and r = p, we express the poten- 
tial as a quartic polynomial function of distance 
and a quadratic function of 



U(r,L, 



GM P 



r 4 
r 4 



B'-± + C-± 



.4 
B 

G 
D 



2Q? 
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To determine the radius of the epicycles (r g ) in- 
duced by a Kepler launch, we foll ow the pro- 
cedure of IS chaffer and Burns (1994), and solve 
for the distance to the potential minimum where 



Or 



= 0. Note that this is only valid 

p—pc ,z—0 

to first order in small quantities, since we are ef- 
fectively assuming that the potential is symmet- 
ric about the equilibrium point. Evaluating the 
derivative, multiplying by r 5 , setting r = + r g , 
and assuming r g << r^, we obtain the epicycle 
radius for a grain launched at in terms of pa- 
rameters known at launch: 



r L (Q, - n L )£l gL 



2 gL -n gL (3Sl + n L ) 



(13) 



In this limit, the radial range of motion of a dust 
grain is simply 2\r g \, and the grain reaches a turn- 
ing point at r t = Tl + 2r g . Note the sign con- 
ventions used here; r g and £l g L may be either 
positive or negative; thus negative grains (with 
QgL < 0) always gyrate to wards R S y n - Equa- 
tion ! 131 corrects a sign error in lSchaffer and Burns! 
(|l994l ) which led to an artificial disagreement be- 
tween the numerical and analytical model in their 
Fig. 6. Equation [T31 by contrast, shows excellent 
agreement with our numerical data for negative 
grains (Fig. [Sk). The peak in Fig. [B^., for oscilla- 
tions towards synchronous orbit, occurs at 



_ r\L / n L -Q 
rg ~ 3 \n L + n 



(») 



n 

II L 



Equation [H] predicts that grains with L* 
launched near R syn reach about halfway to the 
synchronous orbital distance, in agreement with 
Fig. Et. 

For the positive grains, Eq . [TBI gives the proper 
the radial range about stable local minima in both 
the Lorentz limit and in the Kepler limit (Fig.O)). 



9 



At critical values of L*, however, \r g \ —> oo and 
the assumptions under which Eq. [13] was derived 
are violated. This is readily apparent in the de- 
creasing quality of the match between the the- 
ory and the data for intermediate-sized grains in 
Fig. |6b- Note that t his is the same region where 



Mitchell et all (|2003l ) find large non-elliptical gy- 
rations. Nevertheless, the relatively close agree- 



a H (inn) 
0.1 d ^ 




0.01 



Fig. 6. — The radial range of (a) negative and (b) 
positive grains launched azimuthally with the Ke- 
pler speed v — GM/tl at 2.0R P . Both numeri- 
cal data (points) and the analytical results (curves 
from Eq. [T3")) are included. The total radial excur- 
sion is twice the epicyclic radius r g . 

ment between theory and numerical data in Fig. [6] 
confirms that the epicyclic model is usually a good 
assumption in planetary magnetospheres. 

4. Global Radial Stability Analysis 

Our local radial stability analysis makes a num- 
ber of successful predictions, but cannot fully ac- 
count for the boundaries in Fig. [TJd, primarily be- 



cause of the large radial excursions experienced by 
the positive grains. The quartic potential within 
the equatorial plane given by Eq. [T^] contains 
all the information necessary to determine which 
grains strike the planet and which escape into in- 
terplanetary space. 

4.1. Escaping Grains 

Close to the planet, the A/r A term of Eq. [12] 
dominates, and U(r — ¥ 0, L*) — > +oo, while for 
the distant particles we have U(r —> oo,L») — > 
0. Accordingly, the quartic potential can have at 
most three stationary points (one local maximum 
and two local minima). Setting r = gives a 
simple form for the launch potential 



U(r L ,L*) 



GM P 



i 



(15) 



Energetically, a particle is able to escape if 
U(rL,L*) > U(r —> oo,i*) = and we imme- 
diately recover the L* < 1/2 stability criterion of 
lHamilton and Burns! ( 1993bl ). Note that only pos- 
itive grains can escape from a dipolar magnetic 
field and that, in principle, grains with L* > | at 
all launch distances, both inside and outside R syn 
are energetically able to escape. Whether or not 
they do so depends on the form of U(r, L*), in par- 
ticular, on the possible existence of an exterior po- 
tential maximum with U(r pea k, £*) > U(tl, £*)• 

Analysis of Eq.[T2l shows that the potential pre- 
vents all grains launched with Kepler initial con- 
ditions from crossing R syn - Positive grains gy- 
rate away from R syn , while negative grains cannot 
reach R syn (Eq.[f| Fig. [6}. 

Outside R syn , U{r, L») monotonically decreases 
for L* > i. Thus = h is a global stability 
boundary and it matches Fig. QJ> very well. For 
larger L* (smaller grains) , the topography is illus- 
trated in Fig. [7] Stability is determined by the 
height of the distant peak in the potential. For 
L, ~ 1 no such peak exists. For larger L*, how- 
ever, the radial potential decreases with distance 
from tl, then increases to the distant peak, and 
finally declines to zero as r — > oo. 

Consider the quartic equation U(r,L*) — 
U(rL,L*) — 0, which by construction, has one 
root at r = r^, and one root at a more distant 
turning point r — r t . The critical quartic, where 
the turning point is also a local maximum (as in 
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tions where Eq. [12] applies. Since, for positive 
grains in the equatorial plane, the potential de- 
clines as the grain moves inwards from its launch 
distance rt, it can have at most one local max- 
imum within rz,. There are thus two ways in 
which a grain can be prevented from striking the 
planet: i) the potential U at the surface is greater 
than the launch potential, or ii) a potential peak 
exists between the surface and the launch posi- 
tion and its value is greater than or equal to the 
launch potential. These two scenarios are illus- 



Fig. 7. — Potential wells for positive grains 
launched just outside Jupiter's synchronous orbit, 
at r L = 2AR p . For £» = 8.00 {a d = 0.0596jum), 
a distant local maximum bounds the motions. If 
L* = 7.72 (ad = 0.0607^m) the distant peak in 
the potential is at the radial turning point, and 
the potential is equal to the launch potential; this 
is the stability threshold. For smaller L*, the peak 
is lower and escape occurs. 

Fig. [7} has a double root at r = Tf By factor- 
ing out (r — tl), and then differentiating with 
respect to r, we find a quadratic equation for the 
location of the turning point; rt varies smoothly 
from r t 

for r L » R syri 
starts at (r ; 
to 

r L _ 

R.nm. V 27 



?'l at synchronous orbit to r t = |rt 
The stability boundary, tl{L*) 
2 + \/3) and asymptotes 



■R^syn i L 



(16) 



syn • 



syn 

Equation \W\ for r >> R syn is 



for r L » R 
a useful approximation for the boundary far from 
Rsyn, which nicely compliments the exact value 
we have found at the synchronous orbital distance. 
The full solution for the boundary tl(L*) is given 
by a rather messy cubic equation and so we resort 
to numerical methods for its solution, which we 
plot on Fig. [3p. 

4.2. Grains that strike the planet 

Inside R syn , the surface of the planet presents 
a physical boundary to radial motion. Particles 
that strike the atmosphere are slowed and removed 
from orbit. The potential at the planet's sur- 
face U (p, z) varies with latitude, and so for sim- 
plicity, we restrict our attention to planar mo- 



-J 
.j 



-J 




Fig. 8. — Potential wells for two planar trajec- 
tories launched from the solid points which are 
inside Jupiter's synchronous orbit. Distances are 
in planetary radii, and the potential is scaled to 
the launch value. Curve (i) (tl = 1.6i? p ,£* = 
0.0991, a d — 0.53/zm) has a potential peak higher 
than the launch potential inside the planetary sur- 
face. Equating U(rt) = U(R P ) gives an analytic 
solution (Eq. fT7| for the stability boundary in 
Fig.|9t>. Curve (ii) (r L = 1.8i? p ,L* = 0.1277, a d = 
0.48^im) has a potential peak outside the plane- 
tary surface. In this case, the stability boundary 
is best obtained numerically. Both grains depicted 
here are poised on the stability threshold. 

trated in Fig. [SJ For case i), the stability criterion 
is where U(R P ,L*) = C/(r L ,i„). Using Eq.Q^l we 
find a quadratic expression in L* that implies two 
boundaries: 



^-1 
2WRI \R V 



Li 



1 - 



n Lrj 
SIRl 



R„ 



-1=0. 



(17) 
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The two quadratic roots of Eq. [TTl L\ and L 2 , 
may be obtained analytically and are plotted on 
Fig. [9p. The roots obey the simple expression 



(18) 



Equation [18] conveniently highlights several fea- 
tures of the lower curves in Fig.[9p: The two curves 
marking the grains on the threshold of collision 
with the planet are centered on L» < 1, as re- 
quired by Eq. [T5] In addition, for smaller , the 
center of the instability shifts to smaller , hence 
the left-most curve is steeper than the right-most. 
Finally, a planet with a larger R syn (eg. the Earth) 
will have roots that shift to very low L* near the 
planet. 

The curves determined by Eq. [IT] match our 
numerical data cleanly with two important excep- 
tions. Firstly, because our method is only valid for 
grains that collide with the planet in the equatorial 
plane (recall our assumption z = 0), it misses the 
high latitude collisions near [tl — 2R p , L* = 2) 
in Fig. [9]d. All collisions exterior to the bound- 
aries given by Eq. [T7j necessarily involve substan- 
tial vertical motions, and the greyscale shading 
of Fig (9b shows that they do. Secondly, our cri- 
terion predicts instability for a small region near 
(tl = R syn , = 0.2) that our numerical data 
show in fact are stable. These grains encounter a 
high peak, similar to curve (ii) in Fig. that pre- 
vents them from reaching the planetary surface. 
Thus U(rL,L*) > U(R P ,L^) is a necessary condi- 
tion for radial instability in the equator plane, but 
it is not sufficient. 

The additional requirement for instability is 
that U(r L ,L*) > U(r peak , L*), where r peak is the 
location of an interior maximum. Just as for the 
escaping grains exterior to synchronous orbit, eval- 
uation of this condition necessarily involves a cu- 
bic and a semi-analytic method. We find that no 
corrections to Eq. [18] are needed for the high 
radial boundary and for all grains near the planet. 
Only for the right-most curve near R syn is there 
a discrepancy. Our new curve is plotted in Fig [9Jd 
and it perfectly matches the numerical instability 
boundary. Although the stability curve in this re- 
gion can only be obtained semi-analytically, the 
point at which it becomes necessary occurs when 
the potential maximum is located at the planetary 



surface; 



. du 



and U{r L ,L*) = U(R P ,L*). 



r=R„ 



Evaluating these conditions, we find 



(LL. 



I) 5 



,3/2 3/2 



3r L 
Rv 



(19) 



For Jupiter, the critical point that satisifes both 
Eqs. H7] and [H is at L* = 0.112, r L = 1.694R p 
(solid point in Fig[U>). The stability curve meets 
Vl = R syn at L* — 2 — v% a result suggested by 
our local stability analysis of section 3. Note that 
our energy arguments yield analytic expressions 
both inside and outside R syn - Arguments involv- 
ing the location of potential maxima, conversely, 
require semi-analytic methods. 

5. Local Vertical Stability Analysis 

The stability of grains agai nst vertical pertur- 
bation s was first explored by iNorthrop and Hill 
(jl982l ). In their model, a grain is launched on a 
circular orbit at the equilibrium orbital frequency 
u>c in the potential of Eq. @] so that there is no gy- 
romotion around magnetic field lines. If the grain 
orbit at the equilibrium point is stable to vertical 
perturbations, the square of the bounce frequency 
£l b , given by 



n 2 = 



d 2 U 



dz 2 



p—p c ,z=0 



Pi 



9. 



= 3lo 2 c -2u 2 c (20) 



( Northrop and Hill 1982 ) is positive. Here <fi c = 
lj c — Q, is the dust grain's azimuthal frequency in 
the frame rotating with the magnetic field. When 
multiplied by z, Eq.[20l gives the centrifugal (first 
term), and gravitational (last term) accelerations 
along a magnetic field line. 

For < 0, the vertical motion is unstable. 
Note that the gravitational acceleration is nega- 
tive and thus destabilizing. This follows from the 
fact that the dipolar magnetic field curves toward 
the planet and so a grain leaving the equatorial 
along a field line plane moves downhill in the grav - 
itational potential. The lNorthrop and Hill (jl982') 
solution for the boundary where f2& = is plot- 
ted in Figs. QJi and [S^,. At distances closer to the 
planet than a critical distance p cr u , gravity forces 
grains to leave the equatorial plane. 
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5.1. Vertical Instability in the Lorentz 
Limit 



a d (nm) 



0.1 



1.0 



a) 



Lorentz regime 



Kepler regime 




o.oi 



Fig. 9. — Our new analytic results (heavy solid 
lines) are plotted over the numerical data from 
Fig. [T] a) Northrop's solution (dotted line) is 
superseded by our two semianalytic boundaries 
where (f^) = 0% = from Eq. CSD The new 
boundaries are a significantly better fit to the 
data and indicate an inner stability zone. The 
\k c \ = 2f2fc curve indicates the 2:1 resonance be- 
tween the epicyclic and the vertical bounce fre- 
quencies; it matches the data points well, b) Wc 
extend our vertical stability boundary to positive 
grains. The radial stability boundaries for grains 
that escape or crash into the planet are discussed 
in the text (section 4). Between the open circles 
at Tl — R S yn and L* = 2 ± \/3> orbits are locally, 
radially unstable. The solid circle is the critical 
point defined by Eqs. [I~7l and [T9l 



In the limit of high charge-to-mass ratio, Eq. 
can be solved exactly: 



Pcrit 
Rsyn 



(2/3)' 



0.87. 



(21) 



The effect of our initial condition, launching 
grains at the Kepler speed, however, necessar- 
ily causes epicyclic gyromotion as the grain orbits 
the planet. This leads to a stabilizing magnetic 
mirror force, in which the grain resists moving 
out of the equatorial plane to regions of higher 
magnetic field strength as disc ussed in se ction 
1.1. Following the procedure of iLewl ([19611 ) and 
Thomsen and van Alien! ( 1980h , the magnetic mir- 
ror force for equatorial pitch angles near 90° 
adds a component of strength 9r?^ c /2p^ to 
Eq. [20] In the Lorentz limit, Eq. [l3] simplifies 
to rgflg C = p c (fl — n c ), and the bounce frequency 
can be found from: 



ill = 3fi 2 - 2n 2 c + -{n c - ft) 



(22) 



As above, the first two terms are due to the cen- 
trifugal and gravitational forces on a grain tied to 
a nearly vertical magnetic field line. The third 
term of Eq. [22] is the magnetic mirror term, gen- 
eralized to account for a rotating magnetic field. 
The three vertical accelerations add linearly, and 
are valid in the limit that L* — > oo and r g — > 0. 



Fig. [TO] compares the iNorthrop and Hilll (|1982l ) 
bounce period (Eq. |20|) with our Eq. [22] which ac- 
counts for epicyclic motion for small dust grains at 
Jupiter. The Northrop formalism erroneously pre- 
dicts bounce periods that are too long both inside 
and outside synchronous orbit and, more seriously, 
misses the second solution near the planet. 

The third term in Eq. [22] is positive everywhere 
inside the Northrop boundary and thus leads to 
enhanced vertical stability. The stability bound- 
aries in the high \L*\ limit are determined by 
Eq. [25] setting = 0, we find: 



2/3 



0.58,0.84. 



(23) 



These limits are valid for both positive and neg- 
ative grains with \L*\ — > oo. Between these lim- 
its, ill < an d grain orbits are locally unsta- 
ble; the enhanced stability from the mirroring 



13 



3 
C 
J= 




The remaining two accelerations, however, can be 
extended to second order in r g j p c while retaining 
circular gyrations. We begin by writing the verti- 
cal acceleration as a function of the epicyclic phase 



d 2 U _ GM v z{9) ( 3L> 



dz 2 



P 3 (0) 



(24) 



To first order in r„ , the epicycles are circles in the 



guiding center frame. Setting 9 
point to the planet, we find 



at the closest 



Fig. 10. — The bounce period for L* — — 10 4 
grains at Jupiter over a range of launch distances. 
Northrop's solution (Eq. [201 dotted line) and our 
solution (Eq. [251 solid lines) with T b = 2n/U b , 
are plotted alongside numerical data (points). 
Note that in our solution Tf, is smaller than in 
Northrop's solution everywhere except at R syn — 



2.24R p . In the limit where r c << R s 
shows that grains satisfy — > 



n, Eq. HI 

while for 



and Xk — > 3.62 hours. 



force moves the vertical stability boundary in- 
wards from Northrop's 0.87R syn to 0.84i? sjm . A 
more important change, regained stability inside 
0.58 R s yn , is due to the higher launch speeds rel- 
ative to the field lines, larger gyroradii, and a 
stronger magnetic mirror force. For Jupiter these 
distances are at 1.29i? p and 1.87i? p respectively 
(see Fig.[5]). Hints of t his inner stability z o ne were 
seen numerically by Northrop and Hill ( 1983al) 
and iNorthrop and Connernevl (|l987l ); here we 
have derived analytical solutions for vertical sta- 
bility in the Lorentz limit. 



5.2. Vertical instability for all charge-to- 
mass ratios 

To extend our model for bounce motion over 
all charge-to-mass ratios we must, in principle, 
account for the variation in the strengths of the 
vertical gravitational, centrifugal and electromag- 
netic accelerations over one gyrocycle. Extending 
the electromagnetic mirror acceleration requires 
breaking the assumption of perfectly circular gy- 
rocycles, and is beyond the scope of this work. 



p{0) 



and 



Pc 



■ cosO. 



(25) 



(26) 



Due to the geometry of a dipole near its equator, 
an epicycle is tilted by an angle ~ 3A (where A is 
the latitude). Hence the vertical offset is given by: 



z{9) 



3|r 



cos 9, 



(27) 



To calculate the bounce frequency, we average the 
restoring acceleration over an epicycle, a proce- 
dure that is valid as long as n c >> fli,: 



n 2 b = 



, d 2 U 
■ dz 2 



2nz r 



2 7T 



d 2 u 
dz 2 



z(9)d9. (28) 



Using Eqs. [231 and l2"6l to eliminate p(9) and 4>(9) in 
Eq. [Ml we expand to 0(r 2 ), integrate Eq. [251 and 
add in the magnetic mirroring term from Eq. [22] 
to obtain: 



n 2 



3L0 2 -2n 2 c + -(n c -n) 2 



,2 I 2 U 9 



—n, 
2 c 



(29) 



The frequencies in Eq. 



UJ C (Eq.©, n c (Eq.©, 
u gc (Eq. [7]), and 4> c — lo c — f2, are all evaluated at 
the guiding center of motion p c = r^+r g , which is 
determined by Eq. [131 Our calculation adds two 
additional destabilizing terms that are strongest 
for intermediate values of L», where gyroradii are 
largest (Fig. [5]). 

How does our solution compare to numerical 
data? In Fig. we plot our theoretical curves 
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against the numerical data for both negative and 
positive grains launched at the Kepler rate in an 
aligned dipole field for Jupiter. We find the curves 
tracing the unstable zone semi-analytically by set- 
ting fij = in Eq.[25J Within the regions bordered 
by the curves, trajectories are locally unstable but 
may remain globally bound due to high-latitude 
restoring forces. 

Our model closely matches the outer stability 
boundary for negative grains but is less successful 
for the inner boundary, especially for moderate L*. 
This is precisely where our derivation is weakest; 
recall that we have not accounted for higher-order 
corrections for the magnetic mirror force which are 
strongest closest to the planet and for |L*| ~ 1. 
Near = 1 epicycles become large and distorted 
for neg ative grains and even more so for positiv e 
grains (jMendis et all 1 19821: iMitchell etaU l2003h . 
Figure [Bb shows that the epicyclic model matches 
the radial range of positively-charged grains well 
for values L» > 10. This is exactly where the nu- 
merical data depart from the theory in Fig. [9Jd. 
Apparently, large gyroradii and interference from 
the proximate radial instability strip lead to un- 
modeled effects and excess vertical instability. 

The curvature of the outer boundary in Fig. [9£i 
is similar to that for the Northrop instability, al- 
beit displaced to locations closer to the planet. 
Notice that, with decreasing |L*|, the instabil- 
ity region curves towards the planet for nega- 
tive grains, and away from it for positive grains 
(Fig. [9])- This is primarily due to the — 2n^ 
term that determines the Northrop boundary. For 
negative grains inside synchronous orbit, n c > 
uj c , and uj c increases with decreasing |L»| due to 
a weakening outwardly-directed electromagnetic 
force. It thus takes a greater value of n c to make 
3lu c — 2n c change sign, destabilize the vertical 
motion, and move the boundary curves to lower 
launch distances in Fig. |9K For the positive grains 
in the Lorentz limit, by contrast, uj c decreases as 
Li, decreases, and a smaller n c will destabilize the 
grain. Thus with decreasing L* this boundary in 
Fig. curves up to higher launch radii. 

Finally, notice the band of locally unstable but 
globally stable points that stretches from w 
0.1 at the surface of the planet to |L»| w 1 at large 
distances in Figs. QJ, and [9^. These grains are af- 
fected by a |« c | = 2f2b resonance that couples their 
radial and vertical motions. Energy is transferred 



from the radial oscillation to a vertical oscillation 
and back again. Near the synchronous orbit, gyro- 
radii are initially small and therefore there is not 
as much radial motion to transform into vertical 
motion; these grains do not reach our A m = 5° 
threshold and appear as white space in Fig. [9jt. 

The existence of stable trajectories within the 
Northrop boundary is an important result, partic- 
ularly for small slowly-rotating planets with dis- 
tant synchronous orbits like Earth. Small dust 
grains generated by the collisional grinding of par- 
ent bodies on Keplerian orbits can remain in orbits 
near the planetary surface. High energy plasma, 
like that found in Earth's van Allen radiation 
belts, is more stable than we have calculated here 
by virtue of exceedingly rapid gyrations and a 
greatly enhanced mirroring force. 

Our analysis to this point is completely general 
and, although we have focused on Jupiter, can be 
easily applied to other planets. Saturn and Earth 
are logical choices, as their magnetic fields are 
also dominated by the g\o aligned dipolar compo- 
nent. The appearance of the stability map for any 
planet depends on only the parameters R syn and 
R p , and not on the substantially different mag- 
netic field strengths which, due to our use of L*, 
only affect the conversion to grain radius ad- The 
synchronous orbital distance is somewhat closer to 
the planetary surface at Saturn (R syn = 1-&6R P ) 
than at Jupiter (R syn = 2.24-Rp), while at Earth 
(Rsyn — 6. 61 R p ) it is much further away. This 
leads to interesting differences between the plan- 
ets, as we shall see below. 

6. Saturn and Earth 

A centered and aligned dipole is an excellent ap- 
proximation for Saturn's magnetic field. We take 
5io = 0.2154 Gauss from Connernev et all ( 19841 ) 
and plot both our numerical data and analytical 
stability boundaries in Fig. [TT] The Cassini mea- 
surement of <7io does not va ry significantly from 



al. 2009). 



the older value that we use (jBurton et 
A lower synchronous orbit at Saturn pushes the 
local vertical instability inward, as expected from 
Eq. [231 Comparing Fig. [TTJ to Fig. O we see that 
the proximity of the surface at Saturn causes all 
the locally vertically unstable grains to physically 
collide with the planet. This is true for both neg- 
ative and positive grains. 
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Fig. 11. — Stability of charged grains at Saturn 
modeled with a centered and aligned dipole field. 
All initial conditions and theoretical curves are as 
in Fig. IH1 Also as in Fig. the darkest shade of 
grey signifies low latitude collision or escape, the 
middle shade indicates high latitude collisions, and 
the lightest grey signifies large vertical excursions. 
For negative charges (panel a) only a tiny stable 
region exists near (tl = R P , = —50) due to 
Saturn's smaller R syn . Furthermore, due to the 
proximity to Saturn, nearly all grains that are lo- 
cally vertically unstable do in fact hit the planet, 
unlike their counterparts at Jupiter, b) positive 
charges. As with the negative grains, nearly all 
the vertically unstable grains hit Saturn. Saturn's 
radial instability region (darkest grey) looks much 
like Jupiter's. 



Outside synchronous orbit in Fig. [Tib, the solu- 
tions derived for positive escaping grains in section 
4 apply at Saturn to very high accuracy, for both 
the low -L* and high boundaries. As in Fig.O), 
grains with L» < |, do not have enough energy 
to escape despite achieving large radial excursions 
(light grey region outside R syn in Fig. [TTb). For 
these grains, vertical motions are excited over sev- 
eral orbits, as in Fig. 

Within synchronous orbit, the condition U(R P ) - 
U(tl) (solved in Eq. [17)) bounds most of the un- 
stable grains. As at Jupiter, a small set of large 
grains near R syn require the semi-analytical anal- 
ysis of the potential between the launch position 
and the surface to determine global stability. This 
analysis yields the curve connecting the filled black 
circle at (r^ = 1.568i? p , £« = 0.14) and the open 
circle ir L — R syn , L* = 2 — \/3) in Fig. [TTb. 

Compared to Jupiter and Saturn, Earth's mag- 
netic field is "inverted" at the current epoch, with 
magnetic north near the geographic south pole 
(510 = —0.3339 Gauss taken from Roberts and 
Soward (1972). Thus at Earth, L* > for negative 
grains. This causes positive grains to be radially 
stable, gyrating between the launch position and 
synchronous orbit, and negatively-charged grains 
to be radially unstable. The Earth is also far 
smaller on the scale of its own synchronous or- 
bit than the gas giants, and so serves as an ex- 
cellent test of the accuracy of our analytical so- 
lutions far from R syn . For the Earth, Fig. [l"2"k 
shows the radial global instabilities. Outside R syn , 
the boundaries are in excellent agreement with 
our analytical results for large and small grains. 
Inside R syn , grains are globally radially unstable 
and all the grains that collide with the planet at 
low latitudes are launched between our two solu- 
tions given by Eq. 1171 The set of grains for which 
Eq. [T7] is an insufficient criterion for collision with 
the planet, however, is much larger at the Earth 
than at Jupiter or Saturn. For Earth, just like for 
the gas giants, the global solution for radial stabil- 
ity inside synchronous orbit perfectly matches the 
numerical data and meets R syn at the local stabil- 
ity solution (tl = R syn ,L* = 2 — \/3). Further- 
more, the solutions of Eq. [17] have shifted to much 
lower (see Eq. IT8"|). reducing the total range in 
-L* for grains which collide with the planet at low 
latitude. 

The local vertical stability boundary matches 
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the numerical data well, although in the Lorentz 
limit, all grains are globally stable since the high 
latitude restoring forces become much stronger 
close to the planet ([Howard et Only at 

< 1 do the positive grains collide with the 
planet. As in Figs. [9l and ITTI the vertical stability 
curves match very well for large L* and deviate 
from the data for L* » 1. The \k c \ = 2ilf, reso- 
nance also matches the data well. 

Earth has a much larger class of grains that 
experience large radial excursions, which excite 
vertical motions. Most of these grains, from the 
medium-grey areas on the stability map of Fig. 112b 
that link the disjoint dark grey regions of global 
radial instability, collide with the planet at high 
latitudes. An example of a trajectory in this class 
is shown in Fig. [13] 
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Fig. 12. — Stability of charged grains at Earth, 
modelled with a centered and anti-aligned dipole 
field. Theoretical curves and initial conditions 
are the same as in Figs. [9] and QTJ Since R syn 
is much larger than for Jupiter and Saturn, we 
extend the radial range of the integrations to 
tl = 10 R p and the distant threshold signifying 
escape to r esc = 100i? p . The open circles at 
(r L = R syn , L* — 2 + y/3) and (r L = R syn , 
L* = 2 — \/3) are as in Fig. HI and the solid cir- 
cle, marking the transition from the analytical to 
semi-analytical boundary for the larger grains is at 
£* = 0.0248, r L = 2.074i? p . The two solid squares 
in a) arc individual grain trajectories illustrated 
in Figs. [13] and Q31 



Fig. 13. — A grain with large radial excursions 
that gradually excite substantial vertical oscilla- 
tions at the Earth, (r L = 4.51i? p , = 0.948, 
a d = 0.0149/im). 

At Saturn all of the grains in this region col- 
lided with the planet, but at the Earth we see three 
white tracks of orbits that never leave the equato- 
rial plane, and hence are energetically prevented 
from striking the planet. We plot an example in 
Fig. [TJJ A few of these trajectories are also ap- 
parent for Jupiter (Figs. [lb and[3b.) We suspect, 
based on the similarity of the white stable tracks 
in Fig. QUa, and the \k c \ = 2Q b line in Fig. IP2b . 
that these are resonant phenomena. 
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Fig. 14. — A grain with large radial motions like 
that depicted in Fig. Q2] that nevertheless always 
remains near the equatorial plane (rj, = 4.51i? p , 
L* = 0.419, a d = 0.0224^m). 



7. Conclusion 

For Kepler-launched grains in centered and 
aligned dipole planetary fields, we have employed 
both local and global stability analyses to provide 
solutions for stability boundaries that match nu- 
merical simulations for Jupiter, Saturn and the 
Earth. Figure [TBI provides a summary of the var- 
ious analytical results discussed in this work for 
positive grains at Jupiter. 

We find that local radial stability is very use- 
ful in the immediate vicinity of synchronous orbit, 
since r g — > there (Eq. [13]). More importantly, 
our restriction of the global radial analysis to equa- 
torial orbits is justified by the excellent agreement 
between analytics and numerics. Radial instabil- 
ity has important implications for depleting par- 
ticles near the surface of a planet but beyond the 
reach of atmospheric drag forces. At Earth, for ex- 
ample, the radial instability eliminates negatively- 
charged particles with r g < 0.2/^m from Low Earth 
orbit, and < 0.1/xm from within 2000 km. For 
Jupiter, this instability sweeps positive grains with 
r g < 1/im from the region within 10,000 km from 
Jupiter's cloud-tops. 

Our local vertical analysis of grains launched 
on Kepler circles in the equatorial plane adds 
the effect of the magnetic mirror force and is a 
major improvement to t he equilibrium model of 
Northrop and ffill (|l982l ). We do not undertake 
a fully global analysis which would seek to dis- 



Fig. 15. — Local and global stability boundaries 
for positive grains at Jupiter. The dashed line is 
R syn , and the shaded region highlights local radial 
instability near R syn - The white circles bounding 
this region on R syn (L* = 2 ± a/3 for all plan- 
ets) also show where our semi-analytical curves for 
grains that escape or collide with the planet meet 
Rsyn ■ The solid circle is where the potential at the 
planet's surface is a local maximum and equals to 
the launch potential. 



tinguish grains that strike the planet from those 
that simply sustain large amplitude oscillations in 
latitude. 

Although the magnetospheres of Jupiter, Sat- 
urn and the Earth are all nearly dipolar, each 
planet has additional components that make the 
field more complicated. Saturn has the simplest 
field and is well represented by a dipole offset 
northward by a few thousand km. Jupiter and 
the Earth have non-zero dipole tilts that cause the 
magnetic field seen by an orbiting grain to fluctu- 
ate. Nevertheless, since tilts and offsets are gen- 
erally small, we expect that the radial forces will 
be only slightly affected, and the radial instabil- 
ity region will remain nearly the same. Vertical 
motions, by contrast, should be strongly affected 
since a circular orbit in the equatorial plane is 
no longer an equilibrium point. The global radial 
analysis, which included the effects of radial oscil- 
lations, led to a much larger instability region than 
the simple local analysis (top of Fig. fTS"]) ; in exactly 
the same way, we expect the region of vertical in- 
stability to expand substantially when dipole tilts 
or offsets are included. We will take up the study 
of more complicated magnetic field configurations 
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and launch conditions in a forthcoming paper. 
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